Astronomy & Astrophysics manuscript no. 16733 


©ESO 2011 


June 7, 2011 





How to make a clean separation between CMB E and B modes with 

proper foreground masking 

Jaiseung Kim* 



o 

<N 

a 



o 

u 

6 

> 

c3 



> 
m 



Niels Bohr Institute & Discovery Center, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
Preprint online version: June 7, 2011 

ABSTRACT 

We investigate the E/B decomposition of CMB polarization on a masked sky. In real space, operators of E and B mode decomposition 
involve only differentials of CMB polarization. We may, therefore in principle, perform a clean E/B decomposition from incomplete 
sky data. Since it is impractical to apply second derivatives to observation data, we usually rely on spherical harmonic transformation 
and inverse transformation, instead of using real- space operators. In spherical harmonic representation, jump discontinuities in a cut 
sky produces Gibbs phenomenon, unless a spherical harmonic expansion is made up to an infinitely high multipole. By smoothing 
a foreground mask, we may suppress the Gibbs phenomenon effectively in a similar manner to apodization of a foreground mask 
discussed in other works. However, we incur foreground contamination by smoothing a foreground mask, because zero-value pixels 
in the original mask may be rendered non-zero by the smoothing process. In this work, we investigate an optimal foreground mask, 
which ensures proper foreground masking and suppresses Gibbs phenomenon. We apply our method to a simulated map of the pixel 
resolution comparable to the Planck satellite. The simulation shows that the leakage power is lower than unlensed CMB B mode power 
spectrum of tensor- to- scalar ratio r ~ 1 x 10~ 7 . We compare the result with that of the original mask. We find that the leakage power 
is reduced by a factor of 10 6 ~ 10 9 at the cost of a sky fraction 0.07, and that that the enhancement is highest at lowest multipoles. 
We confirm that all the zero-value pixels in the original mask remain zero in our mask. The application of this method to the Planck 
data will improve the detectability of primordial tensor perturbation. 

Key words. (Cosmology:) cosmic background radiation - Methods: data analysis 



1. Introduction 

Over the past few years, CMB polarization has been measured 
by several experiments , and is being measured b y the Planck 
surveyor (iKovac et al.l 
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l2002al; lHa lverson et al.l 
1 2005 1 : lAde & et al.l 120081: 
20091 iBrown & et aP 
Using the general properties 
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Prvke & et al.l 120 09: Hinderks et al.1 
20091: iThe Planck C ollaboration 2006). 

trace-free 



of symmetric trace-free ten- 
sors, we may decompose CMB polarizatio n into the gradient- 
like E mode and the curl- l ike B mode dKamionkowski et all 
119971: IZaldarriaga & Seliakl 1 19971) . In the standard model, B 
mode polarization is not produced by scalar perturbation, 
but solely by tensor perturbation. Therefore, the measure- 
ment of the B mode polarization makes it possible to probe 
\ the u niverse on the energy scale during the inflationary pe- 
riod dKa mionkowski et al .l 119971: [ Zaldarriaga & Seliak] 119971: 
lDodelsonll2003[ lLiddle & Lyth 2000; Mukhanov 20051) . We can 
also obtain information about the dark matter d istribution from 
weak lensing imprints on B mode polarization (lOkamoto & Hul 
120031: ISeljak & Hiratal 12004 iLewis & Challinojl2006l) . In most 
inflationary models, the tensor- to- scalar ratio r is much smaller 
than one, and the current upper bound r < 0.36 i s imposed by 
the WMAP 7 year data at a 9 5% confidence level (lLarson et al.l 
l2010tlKomatsu & et al.ll20l6l) . Therefore, it is quite a challeng- 
ing and ambitious goal to measure CMB B mode polarization. 

Besides contamination inherent to instruments such as 
noise and instrumental polarization, several complications limit 
the detectability of tensor perturbation. In particular, there 
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is contamination from Galactic and extragalactic foregro unds 
dTucci et al] 120051: lAmarie et al] 120051: IVerde et al.1 120061) . By 
spectral matching component separation or template fitting , 
we may clean foreground contamination dBetoule et al] [2009; 
Efstathiou et al. 2009; Bonaldi & Ricciardi 2011). However, we 
need to mask out some regions that cannot be cleaned reli- 
ably. For the analysis of the WMAP data, the WMAP team 
cleaned diffuse foreground by template fitting, and masked 
out 20-30% of a whole sky. Incomplete s ky coverage, due 
to foreground ma sking, leads to E/B mixing(Bu nn et al.l l2003: 
ILewis et aD 120021) . Therefore, there have been various efforts 

to understand and redu ce E/B mi xing a ssociated with in- 
■ * 1 1 — ^ — I I 1 1 1 

complete sky coverage dBunn et al. 1 120031: I Lewis et al. 2002; 
[Lewis 2003; Smith 2006; Kim 2007blal: iKim & Naselskvl feOlO: 
IZhao & Baskaranl 120101: iBunnl 120101) . Previously, we investi- 
gated E/B decomposition in pixel space, and showed that E/B 
mixi ng is localized mostly around the boundary of foreground 
cuts (Kim & Naselskv 2010). 

In real space, E and B mode decomposition operators in- 
volve only the differential of CMB polarization. Therefore, in 
principle, we may succeed in performing a clean E and B de- 
composition in real space. Since it is impractical to apply sec- 
ond derivatives to observation data, we usually perform a spher- 
ical harmonic transformation and an inverse transformation to 
achieve E/B decomposition, instead of applying real-space oper- 
ators. In spherical harmonic representation, jump discontinuities 
in a cut sky produces Gibbs phenomenon, unless spherical har- 
monic expansion consists of infinitely high multipoles. However, 
we may suppress the Gibbs phenomenon effectively by smooth- 
ing a foreground mask, which is similar to the apodization of 
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a foreground mask discussed in other works, though devel- 
oped in a slightly different context ([S mith 200(1 iDas et al.l l2009: 
iKim & Naselskvll20ld: Izhao & Baskaranll2010h . By smoothing 
a foreground mask, we incur foreground contamination, because 
zero-value pixels in the original mask may be rendered non-zero 
by the smoothing process. Therefore, we investigate an opti- 
mal foreground mask, which ensures proper foreground mask- 
ing without any unnecessary loss of sky fraction and suppresses 
Gibbs phenomenon at the same time. We apply our method to a 
simulated map of the pixel resolution comparable to the Planck 
satellite. The simulation shows that leakage power in unmasked 
pixels is comparable to or smaller than the unlensed CMB B 
mode power spectrum of the tensor- to- scalar ratio r ~ 1 x 10" 7 . 
Compared with the result of the original mask, we have reduced 
the leakage power by a factor of 10 6 ~ 10 9 at the cost of only 
a sky fraction 0.07. We find that the enhancement is highest at 
the lowest multipoles and that all the zero-value pixels of the 
original mask remain zero in our mask. 

The outline of this paper is as follows. In Sec.[2l we discuss 
CMB polarization and E/B decomposition. In Sec.[3j we discuss 
E/B decomposition of a masked sky and the Gibbs phenomenon. 
In Sec.|4l we investigate the degree of the Gibbs phenomenon in 
a cut-sky CMB map. In Sec.0 we present a rigorous discussion 
of an optimal foreground mask. In Sec.O we apply our method 
to simulated data and present our results. In Section [71 we sum- 
marize our investigation. 



2. STOKES PARAMETERS 

The state of polarization is descri bed by Stokes parameters 
(lKrausll!986l:lRohlfs & Wilson! 120031) . Since Thompson scatter- 
ing does not generate circular polarization, Stokes par ameter Q 
and U are sufficient to describe CMB polarization (iDodelsonl 
120031) . For the rotation of an angle if/ on the plane perpendicular 
to direction n, Stokes parameter Q and U have the following spin 
±2 properties (IZaldarriaga & Seliak 1997; Zald arriagal 19981) . 



(Q±iUy(n) = e^(Q±iU)(n). 



(1) 



Using its spin properties, we may decompose all-sky 
Stokes parameters in term s of spin ±2 spherical harmonics 
(IZaldarriaga & Seliak|[T997b : 



g(n)±/£/(n) = ^ a ±2 ,i m ±2^/m(n), 



(2) 



/,/// 



where the decomposition coefficients a+2j m are obtained by: 

J [Q(n)±iU(n)] ±2 y; m (n)dn. (3) 



<Z±2,Zm 



Though the quantity shown in Eq.[2]has a direct association with 
physical observables (i.e. Stokes parameters), it is desirable to 
derive rotationally invariant scalar quantities. To derive spin-zero 
quantities , we may apply the spin rais ing and lowering operators 
given by (Zald arriaga & Seljakll 19971) . 
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sin s 6 s mcf>), (4) 
sin'0,/(0,0), (5) 



where s /(0,0) is an arbitrary spin s function. Applying these 
operators to Q(n) ± iU(n) sequentially, we may construct 



the following scalar qua ntities (IKamionkowski et all 1 19971; 
IZaldarriaga & SeljaJdll997h : 



E(n) = 

~[S 2 (2(n) + iU(n)) + 6 2 (g(n) - iU(n))] 9 
B(n) = 

l -[d 2 (£(n) + iU(n)) - S 2 (2(n) - iU(n))]. 



(6) 



(7) 



These two scalar quantities are often termed 'E' and 'B' mode, 
and associated, respectively, with gradient-like and curl-like 
components of the CMB polarization pattern. Besides rotational 
invariance, the construction of scalar B quantities increases the 
detectability of the primordial tensor perturbation, because pri- 
mordi al scalar perturbation make s a null contribution to the B 
mode (Zaldarriaga & Seljak 1997). 

Applying the spin operator s to spin- weighted spherica l har- 
monics, we find the properties (IZaldarriaga & Seljakll 19971) . 

S s Y lm (n) = j(l-s)(l + s+l) s+ iYi m (n), (8) 

S s Y lm (n) = -V(Z + +1) s-iYi m (n). (9) 

Using Eqs.[6l[71[8l and[9j we fin d that the 'E' and 'B' mode quan- 
tities are equivalently given by (Zaldarriaga & Seljak 1997). 



E(n) = 2_j V (i_ 2 )i a EM Y im(n), 
lm y 

B(n) = V J jj^r. a BM Yi m (n), 



where 

^E,lm = 



-(a>2,1m + ^-2,/m)/2, 
K a 2,lm - ^-2,/m)/2. 



(10) 

(11) 



(12) 
(13) 



In real-world observations, angular-scale anisotropy much 
smaller than the FWHM (full width at half maximum) of physi- 
cal beams are sufficiently suppressed. For instance, the window 
function of the FWHM 15 r at the multipole / = 4000 has a value 
~ 10" 24 . Therefore, summations in Eqs.[2l[T0l and [TTl may be 
truncated to a finite multipole with good accuracy. Since it is 
impractical to take second derivatives of observation data, we 
reso rt to Eas.[T0l[TTl[T2l and[T3lto perform the E/ B decomposi- 
tion (IKamionkowski et a l. 1997; Mukhanovll2005l) . 

For a Gaussian seed fluctuation model, the decomposition 
coefficients of E and B mode satisfy the statistical properties as 
follows. 



( a *B,lm aB > l ' m '} = ^l 



EE ( 
BB 



(14) 
(15) 



where (...) denotes an ensemble average. In Fig. [T] we show 
the unlensed Cf E and Cf B values of the WMAP concordance 
ACDM model for various tensor- to- scalar ratio r. As noted from 
Fig. [T] and the WMAP7 upper bounds on r < 0.36, we expect the 
CMB B mode polarization to be much smaller than that of the E 
mode. 



3. E and B mode decomposition of a masked sky 

Several astrophysical emission sources (the "foreground") ex- 
ist between the last scattering surface and our vantage point. 
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Fig. 1. The power spectrum of E and B: B mode power spectrum 
is plotted for various tensor- to- scalar ratio r. 

Contamination from the foregrounds degrades the cosmologi- 
cal information that is attainable from CMB data. Therefore, 
the WMAP team have reduced diffuse foregrounds by template- 
fitting, and blocked some regions by applying a foreground 
mask. That is, reliable measurement of CMB polarization is not 
available over the whole sky, but only a masked sky. We consider 
here E and B mode decomposition from a masked sky. 

E(n) = 
1 



~[S 2 (G(n) + iU(n)) + 6 2 (G(n) - iff(n))], 



B(n) = 

i 

2 

where 



[3 2 (G(n) + iU(n)) - 3 2 (Q(n) - iU(n))], 



(16) 



(17) 



Q(n) ± iU(n) = W(n)(Q(n) ± iU(n)). (18) 

If the sky direction n belongs to the bounded region K of fore- 
ground cuts, a foreground mask W(n) is set to zero, or otherwise 
one. As shown in Eqs. |4] and El a spin raising and lowering op- 
erator are, in fact, second derivatives weighted by trigonometric 
functions. Therefore, a spin raising and lowering operator re- 
quires information within the infinitesimal vicinity of n. Noting 
this, we may show that 

E(n) = -hS 2 W(n)(Q(n) + iU(n))+d 2 (Q(n)-iU(n))], (19) 



[3 2 (G(n) + iU(n)) + 3 2 (Q(n) - iU(n))], 



(20) 



B(n) 



-[3 2 W(n)(Q(n) + iU(n))- 



3 2 (Q(n)-iU(n))], (21) 



= =0 2 (G(n) + iU(n)) - 3 2 (G(n) - iU(n))], 



(22) 



where n belongs to a region outside the foreground cuts and their 
boundaries. In other words, the E and B mode maps constructed 
from a masked sky are identical to those of a whole sky, as far 
as the concerned sky direction does not belong to the foreground 
cuts or their boundaries: 



E(n) = E(n), 
B(n) = B(n), 



(23) 
(24) 



where n £ <R. 

Whether it is a whole- sky coverage or not, it is impracti- 
cal to takes second derivatives of observation data. Therefore, 
we perform a E/B decomposition using Eqs. TUl ITT1 IT2l and [13] 
dKamionkowski et al.|[T997l: lMukhanovll2005l) 7we i 
lently construct E/B decomposed maps by 



may equiva- 



*(») = I J 

lm > 

B(n) = £ J- 

lm y 



+ 2)! 
(1-2)1 

(/ + 2)! 
(l-2)\ 



CiE,lm ^/m(n), 



&B,lm YlmO 1 ), 



where 



&EM - -(ai,lm + 5_2,/ m )/2, 
&B,lm = K^2,lm ~ ^-2,/m)/2, 

and 

a ±2 ,im = J [Q(n)±iU(n)] ±2 7; m (n)dn. 



(25) 
(26) 



(27) 
(28) 



(29) 



As pointed out above, E(n) and B(n) are identical to true E and 
B mode maps for n i Therefore, E and B mode maps con- 
structed by Eqs. [25j [26] [27j and [28] are identical to true E and B 
maps (i.e. those of a whole sky map) for ngR. 



m v V / ( / + 2 ) ! ~ V f N 

E(n) = 2_j ^ (/_2)i a EM y im(v), 

^ /(Z + 2)! 
Bin) = V ( i_2)\ "B,im Y im(ri)' 
lm y 



(30) 
(31) 



4. Spherical harmonic transformation and 
discontinuities 

In deriving Eqs. [30] and [3U we assumed that the forward and 
backward spherical transformations retain a masked polarization 
signal with good accuracy. 



'max 

2(n) ± iU(n) ~ ^ a ±2 ,im ±2 ^/m(n). 

l,m 



(32) 



However, there exists the Gibbs Phenomenon (hearafter GP), 
which refers to the peculiar manner in wh ich the Four i er repre- 
senta tion behaves at a jump discontinuity dGibbsll 1 8981; ICarslawl 
Il92ll) . The Gibbs phenomenon disappear only in the con- 
tinuum limit (i.e. infinite pixel r esolution and summat i on up 
to an infinitely high multipole) (iGibbsl 1 18981: ICarslawl 1 19211: 
lAUen & Mills 2004; Bankman 2000; Proakis & Manolakis 
2006; Glassnei] [l995l) . Jump discontinuities due to foreground 
masking lead to the Gibbs phenomenon, making Eq. [32] a poor 
approximation. To investigate the degree of the Gibbs phe- 
nomenon, we performed a spherical harmonic transformation 
of a masked CMB sky map, and reconstructed it from its in- 
verse transformation. Throughout this paper, all spherical har- 
monic transformation were performed using HEALPix subrou- 
tine map2alm_iterative with three iterations. To construct the 
foreground mask, we combined the WMAP team's polarization 
mask with the point source mask, and prograded it to a higher 
resolution (HEALPix Nside=2048), which is comparable to the 
Planck pixel resolution. The foreground mask described above 
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may be suboptimal for observations other than WMAP (e.g. 
Planck surveyor). However, we find it sufficient for our pur- 
poses, which are the investigation of E/B leakage caused by cut 
sky. For the polarization map, we used the simulated polariza- 
tion map of 10' FWHM, which is described in Section [6j For 
spherical harmonic transformation, we set the maximum multi- 
pole / max to 4096, as recommended by HEALPix for a chosen 
pixel resolution (Nside=2048). In Fig. [2l we show the relative 



separation with proper foreground masking 





Fig. 3. Foreground mask: the input (top) and the reconstruction 
(bottom) 



-0.10 



0.10 



Fig. 2. Relative difference between the input and its reconstruc- 
tion: AQ/Q (top) and AU/U (bottom) 



difference between the input and its reconstruction, respectively, 
for the Q and U signals, where we performed a spherical har- 
monic transformation of the input map and subsequently recon- 
structed it by the inverse transformation. As shown in Fig. [2 
there is a non-negligible discrepancy between the original signal 
and its reconstruction. To show the contours of the foreground 
mask more clearly, we set the minimum and maximum of the 
pixel values to -0. 1 and 0. 1 , respectively, though the value of the 
pixels actually lies in much wider ranges -100 < AQ/Q < 158 
and -160 < AU/U < 502. For comparison, we performed a 
spherical harmonic transformation of Q and U maps without 
the foreground mask and reconstructed Q and U maps by the 
inverse transformation. We find that the relative difference is 
found to be at the level of < 10" 4 . The discrepancy is obvi- 
ously attributed to a foreground mask. For a further investiga- 
tion, we performed a spherical harmonic transformation of the 
foreground mask (HEALPix Nside=2048), and reconstructed 
the foreground mask by the inverse transformation. In a man- 
ner similar to the polarization map reconstruction, we made the 
spherical harmonic transformation up to the multipole 4096. In 
Fig. [3l we show the original mask and its reconstruction. As 
shown in Fig. [3j we find a significant discrepancy between the 
original mask and its reconstruction. In Fig. we show the val- 
ues of the mask W(n) for a fixed azimuthal angle cp = to 
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Fig. 4. W(6, = 0): the red and blue lines correspond to that of 
the input and the reconstruction respectively. 



provide a clearer view, and note the presence of the Gibbs phe- 
nomenon (i.e. the ringing pattern around jump discontinuities). 

5. Smoothed foreground mask and proper 
foreground masking 

The Gibbs phenomenon disappears, when a spherical harmonic 
expansion is made up to an infinitely high multipole. Hence, the 
terms of multipoles higher than the truncation point may be asso- 
ciated with the discrepancy shown in Figs. [2] and [3] Accordingly, 
we may reduce the Gibbs phenomenon by suppressing the terms 
of multipoles higher than the truncation point. In image pro- 
cessing, where the Gibbs phenomenon has been known for a 
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long time, Gaussian smoothin g has been widely used to mit- 
igate the Gibbs phenomenon (|Allen & Mills! 120041; [Bankman 
l2QQQfc iProakis & Manolakisl [20061: lGlassnerlll995l) . To suppress 
the terms of high multipoles, we consider and use the Gaussian 
smoothing kernel consistently throughout this work. 



B(0) = 



1 



V2 



: exp - 



7T<X Z 



2(T 2 9 



(33) 



where 6 is a separation angle. The Gaussian smoothi ng kernel 
is a l ow pass filter, which has the window function (iDodelsonl 
l2QQ3h . 



Bi = exp(-/V 2 ), 



(34) 



where cr = FWHM/ V8 ln(2) and FWHM denotes the full width 
at half maximum of the smoothing kernel. As discussed previ- 
ously, we made a spherical harmonic expansion up to / = 4096, 
as recommended for the HEALPix Nside=2048 pixel resolution. 
We therefore set the FWHM of a smoothing kernel to 15', so 
that the multipoles / > 4096 are sufficiently suppressed (e.g. 
#/=4096 ~ 10" 25 ). One may also consider increasing the maxi- 
mum multipole (i.e. / max ) in a spherical harmonic representation. 
However, given a fixed pixel resolution, increasing / max beyond 
the aliasing limit does not mitigate the Gibbs phenomenon, but 
instead leads to additional distortion. 
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-0.005 0.005 0.01 0.015 0.02 0.025 0.03 

x 

Fig. 5. W(x): the blue, green, red, and cyan curves correspond 
to the original mask, the smoothed original mask, the widened 
mask, and the widened and smoothed mask, respectively. 



Smoothing a foreground mask, we may incur foreground 
contamination, because pixels, which are zero in the original 
mask, may be rendered non-zero by the smoothing process. We 
therefore have to widen the original mask before smoothing so 
that the originally zeroed pixels may remain zero. Without loss 
of generality, we may consider a one dimensional case in a di- 
rection normal to the contours of foreground cut. We assume the 
original foreground mask and the widened mask to be zero, re- 
spectively, for x < and x < xo, as plotted in Fig. [3 where xo is a 
positive number that will be determined later. We apply smooth- 
ing to a widened mask to produce a smoothed mask, which is 
plotted as a cyan curve in Fig. Subsequently, we set the pixels 
(W(n) < P) to zero, where j3 is a very small number. In addi- 
tion to this, we set pixels (W(n) > 1 -J3) to one so that we may 



Fig. 6. Foreground mask: the widened mask (top) and the 'pro- 
cessing mask' (bottom) 

maximize the area of constant values to a good approximation. 
From now on, we are going to call this product 'a processing 
mask' , and use it for E and B mode decomposition via Eq. [2?] 
and [26] To decide the value of J3, we investigated the reconstruc- 
tion accuracy of a smoothed mask, which is mainly limited by 
the numerical accuracy of the HEALPix package. Because the 
associated accuracy is ~ 10" 6 , we set J3 to 10" 6 . As discussed 
previously, we wish to ensure that the zero-value pixels of the 
original mask are zero in the processing mask. To do that, we 
wish the cyan curve in Fig. [5] to have a value of J3 at x = 0, so 
that the pixels x < in the 'processing mask' may be zero. The 
pixel value of a smoothed mask at x = is given by 



= i(l-erf(x /V2^)) 



dx' 



(35) 



where the integral corresponds to the convolution of the widened 
mask with the smoothing kernel. The lower bound of the integral 
is set to xo, since we smooth a widened mask. Using Eq.[35] we 
can show that the boundary location xo of a widened foreground 
mask is given by 



x = V2cr 2 erF A (l -2/3). 



(36) 



However, it is computationally prohibitive to solve Eq.[36l given 
an immense number of pixels. We therefore need to derive 
a computationally efficient form, for which we may consider 
smoothing the original mask. In Fig. the smoothed original 
mask is plotted as a green cyan curve. The pixel value of the 
green curve at xq is then given by 



<2. 



exp - 



Jo 

= ±(l+erf(x /V2^)). 



(x - x'Y 



dx' 



(37) 
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The lower bound to the intergal of Eq.[37]and the exponent argu- 
ment differ from those of Eq.[35l Using Eqs.[36]and[37l we may 
easily show that a - 1-/3, which corresponds to the value of the 
cyan curve at xq. We may therefore produce the widened mask 
efficiently by smoothing the original mask and setting the pixels 
(Win) < 1 - ft) to zero. We emphasize that we have used the 
same smoothing kernel consistently. In other words, cr in Eq.[35l 
should be the same as that in Eq.[37l Implementing the processes 
described above, we produced a widened mask from the original 
mask and a 'processing mask' subsequently from the widened 
mask. In Fig. [6j we show the widened mask and the processing 
mask. We find that the sky fraction of the 'processing mask' cor- 
responds to 0.64. Compared to the original mask (/ S k y = 0.71), 
we find that a decrease in sky fraction is insignificant. We also 
confirm that all zero-value pixels of the original mask are zero 
in the processing mask. 



x 10 
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Fig. 7. W(0, = 0) of the 'processing mask' (a red curve) and 
its reconstruction: The reconstructed 'processing mask' is barely 
visible, because of the visual indistinguishability between the 
'processing mask' and its reconstruction. We also show the orig- 
inal mask (a blue curve), which confirms that the zero- value pix- 
els of the original mask are all masked by the 'processing mask' . 
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Fig. 8. Pixel histogram of difference between the 'processing 
mask' and its reconstruction. 
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AW 

Fig. 9. Pixel histogram of difference between the 'original mask' 
and its reconstruction. 



We performed a spherical harmonic transformation of the 
'processing mask' and reconstructed it by the inverse transfor- 
mation. In Fig. [71 we show the 'processing' mask as a red curve 
and its reconstruction for a fixed azimuthal angle = 0. The 
reconstructed 'processing' mask is barely visible, because the 
'processing mask' and its reconstruction are visually indistin- 
guishable. In Fig. [8l we show the pixel histogram of difference 
between the 'processing mask' and its reconstruction. As shown 
in Figs. [7] and [8j the reconstruction error is quite negligible, 
which indicates that the Gibbs phenomenon is effectively sup- 
pressed. For comparison, we show the pixel histogram of the 
difference between the original mask and its reconstruction in 
Fig.0 



6. Application to simulated data 

Using the WMAP concordance ACDM model, we produced 
simulated Stokes parameter Q and U over the whole sky with 
a HEALPix pixel resolution Nside=2048. For the beam of the 
observation, we assumed a 10' FWHM, which corresponds to 
Planck HFI beam at 100 GHz channel. In Fig.[l0l we show our 
simulated polarization map, where the orientation and length of 
headless arrows indicates polarization angle and amplitude, re- 
spectively. We simulated the inputmap of no B mode polariza- 
tion so that any non-zero values in a decomposed B map may be 
attributed to leakage. For a foreground mask, we use a 'process- 
ing mask' as shown in Fig. We applied the 'processing mask' 
to the simulated map, and produced E and B maps via Eqs. [25j 
[261 [27l and[28l In Fig. [TT] we show the E and B maps, where 
the unmasked pixels correspond to a sky fraction / S ty = 0.64. 
The magnitude of the E and B maps shown in Fig. [TT] are much 
higher than the input polarization map, because of the pref actor 
V(/ + 2)!/(/-2)! in the definitions given by Eqs. [Hand CCD 

From unmasked pixels of the B map, we esti mated the leak- 
age power spectrum us ing the pseudo C\ method (Wa ndelt et al.l 
200 1| iHivon et ai1l2002l) . The power spectrum is usually esti- 
mated using the pseudo C\ method at high multipoles (/ > 30), 
while other metho ds are adopted at the low multipoles (/ < 30) 
(Bo nd et al l 119981: lEriksen et al] 120041: iHinshaw & eTatl 120071: 
Efstathiou 2006). However, we find that using the pseudo Q 
method is good enough for leakage power estimation. In Fig. 
[T2l we show the leakage power spectrum (blue curve) and the 
B mode power spectrum of various values of the tensor- to- scalar 
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Fig. 11. E (top), B map (bottom) 



ratio r. As shown in Fig.[l2l we find that the leakage power (blue 
curve) is comparable to or smaller than B mode power spectrum 
of tensor- to- scalar ratio r ~ 1 x 10" 7 . For comparison, we re- 
peated the analysis, using the original mask (the top figure in 
Fig.O. In Fig.[l2land[T3j we plotted the power spectrum of the 
leakage produced by the original mask as a green curve, and the 
ratio of the blue curve to the green curve. As shown in Figs. fT2l 
and[l3j we reduced the leakage power by a factor between 10 6 
and 10 9 , while losing only a sky fraction 0.07. We also confirm 
that all the zero-value pixels in the original mask are blocked by 
our 'processing mask' . 



7. Discussion 

We have investigated the E and B mode decomposition of 
masked CMB sky data. In real space, the E and B mode de- 
composition operators involve only differentials of CMB po- 
larization. Therefore, we may, in principle, construct E and 
B maps in real space from incomplete sky data without E/B 
confusion. However, it is impractical to apply second deriva- 



Fig. 12. Leakage power spectrum and primordial B mode power 
of various tensor- to- scalar ratio r: blue and green curves de- 
note the leakage power produced, respectively, by a 'processing 
mask' and an original mask. 




Fig. 13. Ratio of the leakage power produced by the original 
mask to that produced by the 'processing mask' 



tives to observation data. Therefore, we usually resort to a 
spherical harmonic transformation and its inverse transforma- 
tion when performing a E/B decomposition. In spherical har- 
monic representation, jump discontinuities in a cut sky produce 
the Gibbs phenomenon, unless spherical harmonic expansion 
consists of infinitely high multipoles. Fortunately, we may sup- 
press the Gibbs phenomenon by smoothing a foreground mask. 



This smoothing approach is simi lar to the apodized foreground 
mask discussed in oth er works dSmithl 120061: |Das et alJ 12009: 
iKim & Naseiskvll2010l: IZhao & Baskaranll201QI) . By smoothing 
a foreground mask, we incur foreground contamination, because 
zero-value pixels in the original mask may be rendered non-zero 
by the smoothing process. We have therefore investigated how to 
derive an optimal foreground mask, which ensures proper fore- 
ground masking without any unnecessary loss of sky fraction 
and suppresses the Gibbs phenomenon at the same time. 

We have applied our method to a simulated map of the pixel 
resolution comparable to the Planck satellite. The simulation 
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shows that the leakage power in unmasked pixels (/ S k y = 0.64) is 
comparable to or smaller than the unlensed CMB B mode power 
spectrum of tensor- to- scalar ratio r ~ 1 x 10" 7 . For comparison, 
we repeated the analysis using the original mask, and found that 
we reduced the leakage power by a factor of 10 6 ~ 10 9 with the 
loss of only a sky fraction 0.07. We also confirm that all the zero- 
value pixels of the original mask are blocked by our processing 
mask. 

Once we construct pseudo E and B maps using the discussed 
method, we will be able to perform analyses such as a power 
spectrum estimation or a Gaussianity study, without worrying 
about the statistical confusion between the E and B mode. The 
application of our work to the Planck data will surely enhance 
the detectability of primordial tensor perturbation. 
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